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We explore the effect of random permanent cross-links on a system of directed polymers confined 
between two planes with their end points free to slide on them. We treat the cross-links as quenched 
disorder and we use a semimicroscopic replica field theory to study the structure and elasticity of 
this system. Upon increasing the cross-link density, we get a continuous gelation transition signaled 
by the emergence of a finite in-plane localization length. The distribution of localization length turns 
out to depend on the height along the preferred direction of the directed polymers. The gelation 
transition also gives rise to a finite in-plane shear modulus which we calculate and turns out to 
be universal, i.e., independent of the energy and length scales of the polymers and the cross-links. 
Using a symmetry argument, we show that cross-links of negligible extent along the preferred axis 
of the directed polymers do not cause any renormalization to the tilt modulus of the uncross-linked 
system. 



PACS numbers: 61.41. +e, 61.43.-j, 62.20.D-, 82.70.Gg 

I. INTRODUCTION 

The statistical mechanics of directed polymers (DPs) 
has been a very active field of research for more than 
twenty years pj [2]. The directed paths under study 
may represent configurations of "real" extended one- 
dimensional objects such as polymers 0] and vortex lines 
in type-II superconductors [I], or may represent config- 
urations in abstract spaces such as those used to model 
sequence alignment in bioinformatics [5j. 

Many physical systems consist of aligned extended one- 
dimensional building blocks which can have crystalline 
or fluid-like order in the transverse plane. Examples in- 
clude columnar phases of DNA [5] , discotic [7J or micellar 
[S] liquid crystals, ferrofiuids 9J, and electrorheological 
fluids [ID] . In addition, polymer brushes consisting of 
dense flexible chains terminally anchored on a surface 
are characterized by chain elongation in the direction of 
the surface normal [TT] . Although the chains of polymer 
brushes can assume backtracking conformations, under 
strong stretching they can be viewed as directed strings 
of Pincus blobs [12 . In recent years, there has been inter- 
est in cross-linked polymer brushes because of promising 
technological applications [121 [H] . 

The effect of quenched disorder in the embedding 
medium on arrays of interacting directed elastic lines has 
led to the prediction of a whole zoo of glassy states in 
high-T c superconductors [4 . In real polymer systems, 
irreversible cross-links can be viewed as quenched disor- 
der of a different type and their effect can be studied 
using the tools of the statistical mechanics of disordered 
systems [15] . A replica field theory has been used to 
study the gelation transition due to permanent random 
cross-links in systems comprised of Gaussian chains |16j . 



beads-and-springs [IT], dimers-and-springs [18] . p-beine 
|19j . and wormlike chains [5DJ. A similar field-theoretic 
approach to well-cross-linked macromolecular networks 
has been developed by Panyukov and Rabin (2TJ [22] 

In this paper, we employ the same theoretical frame- 
work to study the effect of permanent cross-links on a 
melt of flexible directed polymers in a particularly sim- 
ple geometry. The polymers are stretched between two 
parallel flat surfaces with their ends free to slide on them. 
We predict a gelation transition upon cross-linking asso- 
ciated with the emergence of a finite localization length 
in the transverse plane which depends on the distance 
from the boundary surfaces of the slab. 

Furthermore, we investigate mechanical properties of 
the system. Due to the asymmetry of the system, one 
has to distinguish between tilt modulus and shear mod- 
ulus; the first one describes the resistance to shear of 
the boundaries in the preferred chain direction and the 
latter of the perpendicular boundaries. Remarkably, the 
tilt modulus remains completely unaffected upon cross- 
linking with cross- links of negligible extent in the aligning 
direction. 

The paper is organized as follows. We present our 
model in Sec. II. In Sec. Ill, we define and calculate the 
tilt modulus for the cross-linking geometry of our model. 
We construct a replica field theory and obtain the gela- 
tion transition in Sec. IV. The shear modulus is discussed 
in Sec. V. We summarize in Sec. VI. 



II. MODEL 

We consider N directed polymers stretched between 
two planes spaced a distance L apart. The end points 
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FIG. 1: (Color online) A schematic diagram of directed poly- 
mers in a slab of thickness L. z is the preferred direction and 
we refer to x, y as the transverse or in-plane direction. 

of the polymers are free to slide on the planes. Each 
polymer configuration is described by a curve (path) 
r(z) — (x{z), y(z)), where z £ [0,L] and z is the di- 
rection of alignment (Fig. [I]). By the definition of direct- 
edness, these paths exclude loops and overhangs. The 
areal density of the system in the xy plane is N/A. We 
assume free boundary conditions at z — and z = L, 
allowing the polymer ends to assume any arbitrary posi- 
tion on the corresponding planes with any slope. In the 
absence of cross-links, the effective free-energy functional 
( "Hamiltonian" ) of the directed polymers consists of two 
terms, 

+ f 'dzSfaW-r&j), (1) 

where the first term penalizes tilting away from the z 
direction with a being the effective line tension and the 
second term is an excluded volume interaction. 

A physical system where the effective free energy of 
Eq. (JlJ can be realized is that of strongly stretched worm- 
like chains of contour length L. If the two plates are 
held apart by a pressure P and the areal density (on the 
xy plane) of the polymers is high, each polymer will be 
stretched by a tension T ps PA/N. A strong tension al- 
lows only weakly tilting configurations. In this case, the 
aligning part of the free energy reads 




where k is the bending stiffness of the wormlike chains 



related to their persistence length L p via k = L p k B T. 
For T > (k B T) 2 /K and L > L p , one can show [23] that 
the bending term on the rhs of the previous equation can 
be neglected and the projection of the polymer on the xy 
plane behaves as a Gaussian chain. 

In Ref. [12], a realization of a stretched brush is en- 
visioned as ABA triblock lamellae with selective cross- 
linking of the A blocks. In an analogous realization of 
our model, the A blocks would form fluid membranes. 

Another physical realization of Eq. (JlJ is that of worm- 
like chains interacting with a strong nematic field [21 [3] . 
If the chains are embedded in a nematic solvent with 
very large Frank constants and the effective tension 
due to the polymer-nematic interaction a is such that 
a ^> (A;bT) 2 /k, hairpins are negligible. If, in addition, 
L ^> L p , the bending stiffness can be neglected altogether 
and the chains behave as directed flexible polymers. 

In the system described by Eq. |T]), we introduce M 
permanent cross-links which restrict linked polymer seg- 
ments to remain within a distance of order a. Their effect 
is described by an effective interaction 

^=ii>.M-^.)) 2 - <3> 

e— 1 

Cm := {i e ,j e ; z e } is a quenched configuration of M cross- 
links identified by the polymers i e ,j e involved and the 
cross-linking height z e . For the sake of simplicity, we 
assume that the cross-linking interaction is "local" in the 
z direction and only depends on the in-plane distance of 
the polymer segments. As we shall show in Sec. Ill, this 
assumption has profound consequences for the elasticity 
of the cross-linked system. 

The partition function of the system for a specific re- 
alization of cross-links, Cm, reads as 

where (...) denotes average over all polymer configura- 
tions with Boltzmann weight exp(— T-Lo/ksT). Physi- 
cal observables of interest can be calculated from the 
quenched-disorder averaged free energy, F = — k B T [In Z], 
where [...] denotes average over all realizations of random 
cross-links. We assume that the number of cross-links 
can vary and a realization with M cross-links follows the 
Deam-Edwards distribution [T5] : 

The parameter /i 2 = 2[M]/N controls the average num- 
ber of cross-links per polymer, and the physical meaning 
of this distribution is that polymer segments close to each 
other in the un-cross-linked phase have a high probability 
of getting linked. 
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III. TILT MODULUS 



which brings it to the form 



On large length scales, an array of directed polymers 
can be described as an elastic continuum with three elas- 
tic moduli: a shear, a bulk, and a tilt modulus. The first 
two characterize deformations in the transverse plane 
whereas the third characterizes the response to tilting 
away from the preferred axis. The elastic free energy of 
such a system was proposed in the context of vortex-line 
arrays in type-II superconductors by de Gennes and Ma- 
tricon 241: 



d 2 q 



^{(^ + G, 2 )|u(q, fe )| 2 

+ Bq fl q u u^ l (q,q z )u v (-q,-q z )\, (6) 

where K, G, and B are, respectively, the tilt, shear, and 
bulk modulus, fi, v are Cartesian indices in the xy plane. 
u(q, q z ) is the Fourier transform of the two-component 
displacement field which parametrizes the elastic distor- 
tion of the vortex-line array. 

In order to measure the tilt modulus, we consider a 
small force, f, applied at the upper end point of each 
polymer and the opposite force, — f , applied at the lower 
end point. The induced deformation is measured by the 
average tilt field. The tilt field is defined as 



JY 



t(*,*) = £^(* -*(*)) 



(7) 



The energy of the system subject to a cross-link configu- 
ration Cm and to a tilting force f reads as 

Uf=-Ho + H(C M ) - J dz J d 2 r f t(r, z) 
To leading order in f , the average tilt field is 

r-L 



(8) 



1 

~LA 



dz \ d 2 r(t(r,z))i 



K 



f, 



(9) 



where (...)f denotes thermal average with a Boltzmann 
weight corresponding to the energy functional "Hf. The 
tilt modulus, K, can be extracted from the partition 
function as a static linear- response coefficient, 

k n TA 5 2 



1 



where 



K N 2 L Sfntfv 



Z f = J 2?{r 2 (z)}exp^- 



lnZf 



f=o 



(10) 



(11) 



In the path integral of the previous equation, we apply a 
"Galilean" transformation (where the height z is viewed 
as a time-like parameter) [25] . 



r 

st = : 



r 4 + -z 
a 



(12) 



Zf = Z f=0 exp - 



NL 1 

2k B T<r 



(13) 



where Zf = o is the partition function without the external 



field f. Equations (13) and (10) yield 



K 



N 

T 



(14) 



This result implies that the tilt modulus of a directed 
polymer array with a specific realization of cross-links of 
the type described by Eq. ^ is completely unaffected by 
the cross-links and simply reduces to the single-polymer 
tension. Since any realization of the quenched disorder 
associated with the cross-links of this type would give the 
same result, we are spared the burden of having to use 
replicas for the calculation of the average over disorder. 

The reason behind the particularly simple result for 
the tilt modulus is the "Galilean" invariance of the in- 
teractions between the polymers as well as of the bound- 
ary conditions. In the specific model, both the excluded 
volume interaction and the cross-link interaction involve 
polymer segments at the same height z and therefore 
remain unchanged under a "Galilean" transformation. 
Real cross- linking molecules have a finite extent and may 
link to polymer segments at different heights thus break- 
ing the "Galilean" invariance. That would lead to a non- 
trivial renormalization of the tilt modulus. In the limit- 
ing case of cross-links with negligible extent, our model 
is a good approximation, and we expect the tilt modulus 



to remain unchanged and be given by Eq. ( 14 ) 



IV. GELATION TRANSITION 

The system of cross-linked directed polymers under- 
goes a gelation transition as the number of cross-links 
per chain increases. Whereas in the sol phase the DPs 
are free to move in the xy plane like particles in a two- 
dimensional fluid, the polymers' motion in the gel phase 
is restricted to finite excursions around preferred posi- 
tions. Thus there is a localization transition in the xy 
plane, similar to the gelation transition in systems com- 
prised of other building blocks in d = 3 [15] . Since the 
latter has been discussed extensively, we keep our discus- 
sion short. 

What is the order parameter for the localization transi- 
tion in the xy plane? A point z on curve i, i.e., monomer 
z on polymer i in a discretized model, is localized, if it 
has a nontrivial expectation value 



(<5(x-r^)))^l/K 



(15) 



If the particles are localized at random positions, as we 
expect for the gel phase, then the density averaged over 
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all particles vanishes at any nonzero wave vector. A pos- 
sible order parameter is the second moment of the local 
density: 



1 N 

fi(2) (q^)-^E[( eiq ' ri(2) )(^ q - riW )] 



(16) 



In general, one polymer is cross-linked with a finite num- 
ber of other polymers and in fact close to the transition 
this number is small. Hence there is no reason to ex- 
pect that the local density should obey Gaussian statis- 
tics, therefore we need all moments of the local density 
to characterize the gel. This is achieved in the replica 
formalism by introducing n copies, one for each thermal 
expectation value. The order parameter in the replica 
theory 

1 N 

f2(xi...X„,z) = — ^[(d(xi-ri(z))) • • • (6(Xn-Ti(z)))], 



captures all moments of the local density and hence char- 
acterizes the structure completely. 

The average over the quenched realizations of cross- 
links, Cm, is done with help of the replica trick. The 
disorder averaged free energy F — —ksT^nZ] = 
lim n _>o(-2n+i — Z\) / {nZ\) is represented in terms of n 
nonintcracting copies of the system together with one 
additional replica to account for the distribution P(Cm) 
of Eq. ([5]) which is proportional to the partition func- 
tion. The replicated partition function is represented as 
a functional integral over collective fields 0(g, z), 



Z n+ i = / VQe- N ^ , 



(17a) 



■■EL 

2L 



2 pL 



! dzY, \m^)\ 2 ^{q) 



<?ghrs 



l - fdz y, \m^)\ 2 \{q) 



2L .., 

qGIRS 

In 3 , (17b) 



with the single-polymer partition function 



3 = / Vr(z)e~ H o 



(18) 



x exp 



(l f dz S 

\ Jo geUBS 

dz £ A(g)0(g,^)e-^W 



1 

LJ 



where 



(19) 



To simplify the notation we have introduced hatted vec- 
tors, such as q := (qo, qi, ...q„) for (n + l)-fold replicated 
vectors. We have also adopted units of energy such that 
k^T = 1. The harmonic potential for the cross-links is 
reflected in A(q) = exp (-a 2 g 2 /2) and <j) = 2na 2 /A. The 
collective field f2 is almost the order parameter, discussed 
above, except for the zeroth replica which we have intro- 
duced to account for the disorder average [...] in Eq. ( |l~7| . 

Areal density fluctuations are represented by Q(q, z) 
with q = (0, ...q a , ...0), i.e., only one nonzero compo- 
nent (IRS). These fluctuations are penalized by the ex- 
cluded volume interation. The stability of the liquid 
state in mean-field approximation (uniform density) is 
controlled by the coefficient of the quadratic term in the 
fluctuations. A sufficiently strong excluded volume in- 
teraction such that \{q) 3> 1 together with the positive 
definiteness of the kernel (in z\,z%) (e~ l< 7'( r ( Zl ) _r ( 22 ))) = 
exp(— q 2 \zi — z 2 \/2a) preclude a collapse of the liquid 
state. Since the areal density fluctuations are noncriti- 
cal, we only consider the order parameter in the so-called 
higher replica sector (HRS) consisting of vectors q with 
at least two nonzero components. 

The expectation value of the order-parameter field 

(Q(x,z)) f = (6(x-f(z))) f (20) 
has to be calculated self-consistently with the weight of 



Eq. (17). Here we restrict ourselves to the saddle-point 



approximation 8f n +i/Sfl — 0. As for the gel transition of 
random coils, the saddle-point equation is solved exactly 
by the following ansatz for the order parameter: 



il(q,z) = (l-Q)8 9fi 

+ QS^, J dli 2 V(e,z)exp(- 



(21) 



Here Q denotes the fraction of DPs in the infinite clus- 
ter and hence 1 — Q is the fraction of DPs in the fluid 
state, giving rise to the first (trivial) contribution to the 
order parameter. On the other hand, the localized parti- 
cles are characterized by the localization length £, which 
fluctuates not only from polymer to polymer but also 
along one directed polymer giving rise to a distribution 
of localization length 



QN 



(22) 



which depends on z. After averaging over the disorder 
the system still has macroscopic translational invariancc 
in the xy plane. This requires that := J2'a=o I" = 0- 
The solution of the saddle-point equation reveals a gel 
transition at a critical cross-link concentration fi = 1, 
when a macroscopic cluster of cross-linked DPs is formed. 
The gel transition is signaled by a nonzero value of the 
gel fraction Q. Close to the gel transition, f = (i 2 -l«l 
grows continuously from zero: Q — 2e+C(e 2 ). To discuss 
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the distribution of localization length, we note that we 
have several length scales in our system: the internal 
length of a directed polymer, L, the length of a cross-link, 
a, and the radius of gyration in the xy plane, which is 
determined by I :— y/L/ (2a). We expect that the latter 
will set the scale for the localization length, introduce 
the abbreviation 9 = (2/3 + a 2 / l 2 )l 2 / (t; 2 e) and consider 
the distribution of rescaled, inverse localization length, 
ir(9, s — z/L), with 



n(6,s)de = V(£, 2 ,Ls)d£, 2 
This function is the solution of 

2e)7r(0, s) = (1 + e) / d Sl ir(9, Sl ) 



(1 



2/3 + a 2 /P 



ds 1 d e (e\(e,si)){2\ t 



si 



(23) 



(24) 



7H 



dsids2 / d9\ 7r(#i, Si)ir( 



s 2 )+0(e 2 



To gain a better understanding of the solution, we 
decompose the distribution into its mean with respect 
to s, 7f(#) = J ds 7r(#, s) and a deviation: 7r(0, s) = 
n(9) + dir(9,s). The mean, n(9), fulfills the same equa- 
tion as for isotropic gels [TB] . The deviation is small close 
to the gel point, 



6ir(9,s) = ew{s)de(9 2 Ti(9)^ + 0(e 2 ) , 
s 2 + (1 - s) 2 - 2/3 



w(s) 



a 2 /I 2 + 2/3 



(25) 
(26) 



and furthermore controlled by the ratio of cross-link 



length to in-plane radius of gyration: 



a 2 /I 2 . The 



larger the radius of gyration I, the more pronounced is the 
dependence on s. We show the distribution for a typical 
value a 2 = 0.1 in Fig. [2] As one would expect, localiza- 
tion is strongest in the middle of the directed polymer 
and weaker at the boundaries. The variation across the 
length of the DPs is stronger for larger localization length 
(small 9). 

To get a better understanding of this anisotropy, we 
show in Fig. [3] cuts of Fig. [2] for two fixed 9 values. As 
one can see, large localization lengths (such as 9 = 1/2, 
solid red curve) are favored at the boundaries (s 0, 1), 
and small localization lengths (such as 9 — 2, dashed 
blue curve) in the middle of the sample (s sa 1/2). This 
behavior is reasonable, since the ends of the chains are 
more loose. A chain segment close to the top (bottom) 
boundary has a lower probability to have a cross-link 
above (below) and hence is on average less localized than 
a chain segment in the middle. 



V. SHEAR MODULUS 

In the gel phase, the DPs are localized in the xy plane. 
Hence the symmetry with respect to translations in the 
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FIG. 2: (Color online) Height-dependent part, 5tt(8, s), of the 
distribution of inverse localization length, 9, as a function of 
9 and s := z/L. Due to the scaling behavior close to the sol- 
gel transition, 5ir(8,s) is normalized with e. af is the ratio 
squared of the in-plane extent of a cross-link to the in-plane 
radius of gyration of a free DP. 
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FIG. 3: (Color online) Height-dependent part, 5-k(9, s), of the 
distribution of inverse localization length, 9 vs. the normalized 
system height s — z/L, for 9 = 1/2 (solid, red) and 9 = 
2 (dashed, blue). The cross-link extent is a 2 — for both 
graphs. As in Fig. [2] 8tt(9, s) is scaled with e. 



x and y directions is spontaneously broken as indicated 
by a nontrivial expectation value of the local density as 
defined in Eq. (15 1. The symmetry breaking occurs on 



a local level only, while the macroscopic system (aver- 
aged over the cross-link related disorder) remains homo- 
geneous. We expect low-energy Goldstone fluctuations 
and a finite stiffness to static shear deformations in the 
xy plane (see Fig.[T]). In the replica formalism, the overall 
macroscopic translational invariance is reflected in com- 
mon translations of all n + 1 replicas, whereas replica- 
dependent translations generate a family of order param- 
eters which all give rise to the same free energy. 

To investigate the response of the system to shear 
deformations, we start from the replica free energy 



[Eq. (17b)] and consider fluctuations around the saddle 
point (21), which correspond to long- wavelength shear 
deformations u"(r) in each replica a = l,...,n except 
the zeroth replica [T71 [27] . The latter represents the 
preparation state before the cross-linking process, which 



(> 



takes place in a state without shear deformations. We 
furthermore want to consider pure shear only and hence 
require that volume is conserved, V ■ u a = 0, and no 
tilt deformations are excited; i.e., u"(r) are chosen to be 
independent of z. 

The order parameter for the deformed state thus reads 

as 

Q u (q,z) = (1-Q)S~ Q + 

-^exp(iq r r + ^qi-u«(r)) 



a=l 

ee 



(27) 



where qy"' : 



,(«) 



If the deformations are taken 



to be spatially uniform we recover the general solution 
of the saddle-point equation. Fluctuations around the 
saddle-point value are incorporated by nonzero d x u and 
d y u. These are assumed to be small, corresponding to 
long-wavelength excitations. 



We plug the ansatz (27 1 into the free energy ( 17b) and 



only keep the lowest order in Q and in the derivatives 
d x u and d y u. Higher order derivatives, such as d^ y u, 
are neglected. The result has the form of an elastic free 
energy of an incompressible medium, 



/n+i(0„) = U + ^ fd 2 rJ2Y, 

J , , , 1 „ i 



dr u 



(28) 



There / sp is the saddle-point value of the free energy and 
G is the shear modulus, 



2 6 J A 



(29) 



With the distance from the sol-gel transition e = /x 2 — 1, 
the relation Q — 2e + (D(e 2 ) found in the previous section, 
and the areal density no := N/A of the polymer chains, 
the shear modulus simplifies to 



G = -e 3 n Q k B T 



(30) 



The scaling of the shear modulus Gat 3 close to the 
sol-gel transition is in agreement with previous results 
for isotropic systems [T71 [551 |S7]. This result for the 
shear modulus is universal and does not depend on the 
microscopic length or energy scales which characterize 
the polymers and the cross- links. 



VI. CONCLUSIONS AND OUTLOOK 

We have addressed the effect of random permanent 
cross-links on an array of directed polymers confined be- 
tween two planes with their end points free to slide on 



them. The cross-links are assumed to have negligible ex- 
tent in the preferred direction of the DPs (local in the 
z direction), but they are spring-like in the transverse 
(xy) plane. The constraints imposed by the cross-links 
are treated as quenched disorder which follows the Deam- 
Edwards distribution. 



At a certain critical cross-link density, there is a con- 
tinuous gelation transition from a sol phase where the 
DPs are free to wander in the x and y directions to a 
gel characterized by the emergence of finite localization 
lengths for the polymer segments which belong to the in- 
finite percolating cluster. Unlike other isotropic polymer 
systems which undergo a similar transition, the DPs are 
inherently anisotropic and this is reflected in the height 
(z) dependence of the order parameter and the associated 
distribution of localization length. Because of the finite 
extent of the system in the preferred direction, larger lo- 
calization lengths are favored closer to the boundaries 
where the polymer end points are free to slide. 



The gelation transition is accompanied by the emer- 
gence of a finite shear modulus. Our result for the array 
of cross-linked DPs close to the gel point agrees with pre- 
vious results for isotropic systems thus suggesting uni- 
versality. As far as in-plane localization and the rele- 
vant shear modulus are concerned, our system can be 
viewed as effectively two dimensional. It is well known 
that truly long-ranged positional order cannot exist in 
two dimensions [2 . In [27 , it was shown for isotropic 
systems that fluctuations drive the order parameter to 
zero as expected from the Mermin- Wagner theorem. Yet 
a quasi- amorphous solid state survives. It is character- 
ized by a finite stiffness to static shear deformations and 
algebraically decaying correlations. 



The asymmetry of our system due to the preferred di- 
rection of the DPs entails the existence of a tilt mod- 
ulus which is different and independent from the shear 
modulus. We have only considered the simplest case of 
cross-links with negligible extent in the preferred direc- 
tion of the DPs. Using a "Galilean" invariance argu- 
ment, we have shown that cross-links of this type leave 
the tilt modulus of the uncross-linked system completely 
unaffected. We expect cross-links which connect poly- 
mer segments at different heights to induce an effective 
"nonlocal in z" interaction between the connected poly- 
mers. By analogy to a similar interaction in the case of 
flux lines in type-II superconductors [2S], we can expect 
the breaking of the "Galilean" invariance to cause an 
upward renormalization (stiffening) of the tilt modulus. 
This putative renormalization may be useful to quantify 
the cross-link induced collapse of polymer brushes. We 
hope to report on these issues in a future publication. 
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